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Analytical theory of fast magnetic reconnection with a large guide field is presented for 
the first time. We confirm that two distinct reconnection regimes are possible depending 
on whether the diffusion region thickness S is larger or smaller than the sound gyroradius 
p s . The reconnection is slow or Sweet-Parker-like for 8 > p s , and fast otherwise. In the fast 
regime, however, we find that ion viscosity p plays a critical role. In particular, for 5 < p s the 
diffusion region thickness is proportional to Ha^ 1 with Ha oc \j ^fffp the Hartmann number, 
and the reconnection rate is proportional to Pr -1 ' 2 with Pr = p/n the Prandtl number and 
rj the resistivity. If the perpendicular ion viscosity is employed for p the reconnection rate 
becomes independent of both plasma j3 and collision frequencies and therefore potentially 
fast. 



Fast magnetic reconnection phenomena observed in space and laboratory plasmas indicate that 
nonlinear reconnection rates are independent of collisional dissipation coefficients. Such dissipation 
independent rates have been reproduced computationally for reconnection with lH5fl and without 
Q| a guide magnetic field. They have also been predicted analytically in electron 7| and Hall 10 1 
magnetohydrodynamics (MHD) models without a guide field. However, no analytical model for 
reconnection with a guide field is currently available. Since a guide field is normally present and 
sometimes significantly exceeds the reconnecting in-plane magnetic field components, as in magnetic 
fusion experiments, a nonlinear analytical model of fast reconnection with a guide field is highly 
desirable. 

Extended MHD simulations of guide-field reconnection have demonstrated that fast reconnection 
regimes [l| 



can be enabled by the electron pressure gradient in the Ohm's law [ll| and occur 
when the sound gyroradius p s , defined as the ion gyroradius with electron temperature, exceeds the 
diffusion region thickness 



Ea|. 



[3Q. 



12[. The latter conclusion was recently confirmed experimentally 



However, despite significant theoretical progress, a number of fundamental questions remains 
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unanswered. In particular, the scaling of the nonlinear reconnection rate with p s is not unique and 
appears problem dependent: Kleva et al. [1] found E z oc p® for magnetic flux bundle coalescence; 



n 



Bhattacharjee et al. |l4j | observed E z oc p s for tearing modes; Fitzpatrick Q| concluded that 
E z oc p'J' 1 for the Taylor problem; while Schmidt et al. [5] saw E z oc /?" for tearing modes and 
magnetic island coalescence, where a decreased from 1 to as p s increased. Formation of a thin 
current sub-layer inside the current sheet layer, which results in quasi-explosive reconnection during 
early nonlinear stages 



14 



18], and its nonlinear saturation are also poorly understood. Simulations 
show that the total layer thickness remains of order the electron skin depth d e [lfil. \\\. whereas the 

16l4l8| (but not demonstrated) 



sub-layer thins with time without a limit. It has been speculated 
that such thinning in unphysical and will eventually be arrested by three-dimensional nonlinearities 
or electron viscosity p e . Ion viscosity p is usually believed to play no significant role in the large- 
guide- field fast reconnection process (e.g., see Refs. flfiEJ). 

Here we propose, for the first time, a simple nonlinear analytical model of magnetic reconnection 
in extended MHD with a large guide field. We consider a quasi-steady state diffusion region in two 
dimensions and evaluate its aspect ratio and the corresponding reconnection rate. We take into 
account plasma resistivity 77, ion viscosity p, as well as finite sound gyroradius p s effects. Similar 



to the zero-guide-field case [7Bl0j|. we find that two reconnection regimes are possible: a fast and 
a slow, Sweet-Parker- like [20]. But, unlike the zero-guide-field case, where the reconnection rate in 
the fast regime is not explicitly dependent on dissipation coefficients, here we find it proportional 



to Pr 1 / 2 with Pr = p/r] the magnetic Prandtl number. Assuming 2lJ that p corresponds to 



the perpendicular ion viscosity 



22| gives Pr" 1 / 2 ~ /3- 1 / 2 (m e /m i ) 1 / 4 (r i /T e ) 1 / 4 , which predicts the 



reconnection rate not to be explicitly dependent on collision frequencies. Here /3 <C 1 is the ratio 
of plasma pressure to magnetic pressure, and m e , nii, T e , and Tj are electron and ion masses and 
temperatures, respectively. Therefore, contrary to previous beliefs [l[ ^,0], ion viscosity plays a 
fundamental role in the fast reconnection process with a guide field. A transition between the two 
regimes occurs when the current layer thickness becomes comparable with p s , in agreement with 
previous findings [l|, [5], Q] • The model gives an explicit dependence of E z on p s and suggests 
why different scalings can be observed for different reconnecting systems. It also obtains a nonlinear 



18| and demonstrates how it is arrested by dissipation 



threshold for the current sheet thinning [14 
processes. 

Two-field fluid model. - We consider plasma in a strong, straight, homogeneous magnetic field 
Bq = BqZ, z = Vz, and employ the two-dimensional (d/dz = 0), low-/3, cold-ion fluid model of 



Refs 



14, 



17 



21 



, [23j] for magnetic reconnection in the x — y plane. For simplicity, we neglect in 



this work finite electron inertia d e effects. These effects are important and will be discussed in future 



work. Normalizing to the Alfven speed Va = Bq/ V '4-Knomi, with no the homogeneous equilibrium 
plasma density, and an arbitrary length L gives 

d 



at 



+ V • V \vj = B ■ V(V-ty) + pV 2 w, 



dB 

~dt 



+ V x [B x (V - p\z x Vw)] = -7]V x (V x B). 



(1) 
(2) 



Here, V = z x V(/? is the ion flow velocity with ip the electrostatic potential, w = V tp is the 
flow vorticity, and £? = i: x VV' is the magnetic field in the x — y plane. The sound gyroradius is 
p s = \jT e jmil '(£liL) with 0^ = ei?o/(^jc) the ion gyrofrequency, c the speed of light, and e the 
proton electric charge. Finally, rj and p are normalized resistivity and ion viscosity, respectively. In 
numerical simulations, Eq. ([2|) is often replaced with a single scalar equation for ip. However, for 
our purposes [^-[sl. 



it is convenient to consider the vector form of the equation. 




► W 



Figure 1: (color online). Diffusion region geometry. 



We concentrate on the diffusion region and do not concern ourselves with a particular physical 
system supplying magnetic flux. A closed dynamical description can be obtained b y c oupling our 



microscopic diffusion region description with a suitable macroscopic driver, as in Ref. [24j |. As shown 
in Fig. Q] we assume a rectangular diffusion region of (normalized) dimensions 5 and w with plasma 
entering and exiting along the y = Vy and x = Vx directions, respectively. We define the discrete 
upstream and downstream magnetic field variables B x = x ■ B(0,5/2) and B y = y ■ B(w/2,0), 
respectively, and the discrete flow stream function <I> = —ip(w/2, 5/2). Then, the inflow and outflow 

. nun 

velocities are given by V y = —2&/w and V x = 2<3?/<5, respectively. Following Refs. [7-9, 19] we next 
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discretize Eq. ([T]) at (x,y) = (w/2,5/2) and the x and y components of Eq. ([2]) at (x,y) = (0,5/2) 
and (w/2,0), respectively. Time-derivative terms are normally found small in nonlinearly saturated 
states, e.g. at and around the time of the local maximum of the reconnection rate (2^. A steady- 
state analysis of the diffusion region (but not the full domain) is appropriate in such situations. 
Neglecting time derivatives and numerical factors of order unity and following a similar procedure 
as in Refs. [tHI], results in the following equations for <3?, B x , and B y : 

£(^-i) + (^ + T)(^- ; r)=-'" ,, (i + i) 2 ' (3) 



Sw 
$B 



1 + ^(7^ + 



1 | 1 
S 2 w 2 



!/ 



5w 



1 + Ps(t^ + 



1 1 



5 2 w 2 



Sw 5 2 ' ' 

ow w z 



Equations ([3]) - ([5]) are invariant under theplasma flow reversal, i.e., under the substitution 
(B x , By, 5) (By, — u>), as expected [7|,ll9||. These three equations contain five unknowns: 
5, w, B x , By, <£. It is therefore necessary to consider two independent parameters. We have found it 
convenient to choose w and B x as parameters, as they can be found via coupling with a macroscopic 
driver Q. 

Solution of the discrete equations. - An important characteristic of the diffusion region is its 
aspect ratio £ = S/w. The reconnection rate E z = —r\J z = n(B x /5 — B y /w) can be conveniently 
expressed in terms of £ as 

_ e z _ V2(i-e) ftt . 

z * - B 2 ~ S V C ' 1 ' 

where S„ = \f2B x wjr\ is the resistive Lundquist number. It is clear from Eq. §G§ that, for given 
B x and w, large reconnection rates preferentially occur for £ <C 1. We will concentrate next on the 
£ < 1 limit by approximating 1 + £ 2 ~ 1 — £ 2 ~ 1 in Eq. ^ and elsewhere. Then, the reconnection 
rate becomes E z * ~ y / 2/(5' r? £). An equation for £, or equivalently 5, can be obtained from Eqs. ([3]) 
- ©. Introducing a dimensionless quantity 5 = S/p s , which characterizes importance of plasma 
compressibility, and eliminating variables results in the following equation for £: 

1 , 1+^ _ (i+jr (7) 



where S 1 ^ = \/2B x w/fi is the viscous Lundquist number. 

To solve Eq. ([7]) for £, or equivalently 5 = £w, it is convenient to introduce parameters q = 
Pr/(Pr + 1) < 1, with Pr = fi/rj the magnetic Prandtl number; and 5sp = Ssp/Ps, with 5$P = 



5 



(w I \fS^){\ + Pr) 1 / 4 the Sweet-Parker length scale. Then, Eq. ([7]) becomes 



5 2 (5 2 + 1) 



5 i SP (5 2 + q) = 0, 



so that 



Ssp, $SP > Ps 

k VQ^sp/Ps = w 2 /(p s Ha), 5 S p < p s 
where .ffa = ^S^S^ is the Hartmann number. The corresponding reconnection rates are 

V2/S^(Pr + I)" 1 / 4 , <5 5 p » p s 



(8) 



(9) 



(10) 



Consequently, two reconnection regimes are possible. When the diffusion region thickness exceeds 
p s the reconnection is Sweet-Parker like and therefore slow. In the opposite limit the diffusion 
region thickness is proportional to Ha -1 and the reconnection rate to Pr -1 / 2 = \Jt)I \i. Therefore, 
both ion viscosity and resistivity are essential ingredients. However 2l|. if p corresponds to the 
perpendicular ion viscosity, then Pr~ 1 / 2 ~ /3~ 1 / 2 (m e /mj) 1 / 4 (Tj/r e ) 1//4 and the reconnection rate is 
not explicitly dependent on collision frequencies. Hence, it is potentially fast. Since p s = yf3/2di 



with di = (c/ ' L) nii / (4-Kne 2 ) the normalized ion inertial length scale, the reconnection rate E Zif = 
(2/Pr) 1 ^ 2 (p q /w) in the fast reconnection regime is independent of <C 1. A transition between the 

nnnn 

two regimes occurs at 5 ~ 5sp ~ p s , as expected [1|, |5|, H2|, ll3fl. These conclusions are confirmed 
by solving Eq. ([8]) numerically. The solution for q = 0.01 is shown in Fig. [2] with a solid line. The 
asymptotics ([9]) are shown with dashed lines. 




Figure 2: (color online). Numerical solution of Eq. JHJ with q — 0.01 (solid line). Asymptotics ((9]) are shown 
with dashed lines. 
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It follows from Eqs. ([9]) and (|10p that the fast reconnection regime does not exist in steady state 
for p = 0. This corresponds to the absence of a steady-state solution and the indefinite current 
sub-layer thinning observed in Refs. |l4 - 



18| . In fact, Eq. ([8]) predicts that a steady-state solution 
in this case only exists for 6sp > p s - Finite ion viscosity regularizes the solution for 5,dsp < P- 
Similarly, Eqs. ([9]) and (jlOp predict that the fast reconnection regime does not exist when p s = 
for arbitrary p, as expected. 

Fundamental role of ion viscosity. - The relevance of p in the fast reconnection regime with a 
large guide field can be understood directly from Eqs. ([1]) and ([2]). Consider a quasi-steady state 
diffusion layer and assume £ <C 1, so that according to Eqs. ([4]) and ([5]) B y /B x = (< 1. Equation 
dU) requires the term B • V(V 2 ^) to be balanced by either V ■ Vro or pSJ 2 w. Therefore, viscosity 
is important if pV 2 ~ p/8 2 > V • V ~ $>/(5w) and unimportant otherwise. The magnitude of $ 
can be estimated from Eq. ([2]) by discretizing its x component at (x,y) = (0,5/2): 

B x <& ( p 2 \ riB x 



<5u> V <5 2 / 5 2 
Then, 

• ~ I "' S2 - S >P - . (12) 
^ 1 Til . S < ft 

Since rj and /i have the same scaling with collisionality when p represents perpendicular ion viscosity, 
and since 5 decreases with collisionality, for 5 < p s the viscous term in the vorticity equation will 
dominate over the advective term for small collision frequencies and p will play a fundamental role. 
On the other hand, for 5 > p s both terms are comparable and p does not influence reconnection 
fundamentally. Employing result (|12p in Eq. (JTJ) discretized at (x,y) = (5/2, w/2), 

$ 2 V$ Bl 



5 3 w + 5 4 ~ Sw' (13) 

gives result ([9]) for 5 and consequently result (jlOp for E z *, as expected. 

Comparison with previous numerical studies. - The results obtained herein agree well with the 
available numerical simulations of guide field reconnection. Specifically, Kleva et al. [1] employed 
Eqs. (H|) and ([2]) with p = r] (i.e., Pr = 1) to study coalescence of two magnetic flux bundles. The 
work varied p s and rj in the fast reconnection regime (<5 < p s ) by as much as a factor of four each 
and found w oc p s and E z oc rj ,p®, consistent with Eq. (|10p . In addition, it was observed that 
5 oc 77, in agreement with Eq. ([9]) when p = rj. 
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Schmidt et al. Q| also employed Eqs. ([1]) and ([2]) with p = r] (again, Pr = 1), but studied recon- 
nection driven by a tearing mode and by a magnetic island coalescence instability. In both cases, 
the reconnection rate in the fast reconnection regime (5 < p s ) was again found to be independent 
of r/, E z oc rf . Moreover, it was observed that E z oc p" with a ~ 1 for 5 ~ Ssp ~ p s (i- e -j at the 
slow-to-fast transition) and a ~ for p s 3> 5s p- Both scalings are also consistent with Eq. (1 1 [) 
when one considers that, at the transition, w is still determined by the Sweet-Parker dynamics and 
is therefore independent of p s (i.e., w oc p®; hence E z oc p s /w oc p s ), while w oc p s when 5 is deep 
in the p s sub-layer [l|, resulting in E z oc p s /w oc p®. 

One set of results apparently at odds with our predictions is documented in Ref. [3], where a 
system of four nonlinear equations for tp, tu, and the z components of the perturbed magnetic field 
and ion flow velocity is employed that is valid for both [3 <C 1 and /3 3> 1 regimes. When solving this 
system in the nonlinear reconnection regime with a large guide field, it was found that E z oc rp p^J 2 
with fixed p, which appears to disagree with our findings ([9]) and ()10p . However, no data is provided 
for w in the reference. Furthermore, these four equations can only be reduced to our Eqs. ([1]) and 
((2|) when [3 <C 1, k = p <C f3ru and v = 0, with k the plasma heat conductivity and v the hyper- 
resistivity as defined in Ref. Since Ref. Q| employed k = p ~ fir] and v = 2.5 x 10~ 9 > 0, the 
results therein cannot yet be compared with our theory. 

In conclusion, we have presented a simple nonlinear analytical model of magnetic reconnection 
with a large guide field. We found that two reconnection regimes are possible depending on whether 
the diffusion region (current layer) thickness 5 is larger or smaller than the sound gyroradius p s . For 
S > p s the reconnection is slow and is described by the standard Sweet-Parker expressions modified 
to account for ion viscosity p. For 5 < p s the character of reconnection changes. Instead of scaling 
with S-q the diffusion region thickness becomes proportional to Ha~ l with Ha = ^JS^S^ the 
Hartmann number; and the reconnection rate becomes proportional to Pr -1 / 2 with Pr = pjr\ the 
Prandtl number. Assuming that p describes the perpendicular ion viscosity results in a reconnection 
rate that is independent of collision frequencies and plasma j3 and is therefore potentially fast. The 
transition between the two regimes occurs at 6 ~ Ssp ~ p s - Therefore, contrary to the common 
belief , ion viscosity is found to play a fundamental role in large-guide-field reconnection. In 

particular, a steady-state fast branch does not exist for p = 0, resulting in indefinite thinning of 5 
with time fl^iO] when 5 < p s . 

Given the demonstrated relevance of ion viscosity for the low-/3 Hall MHD reconnection, it seems 
clear that the simple resistive and viscous closures commonly used in the literature for Eqs. ([TJ 
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and ([2]) are likely to be inadequate for accurately describing weakly collisional magnetospheric 

ows that more accurate models for collisionless electron and 



and magnetic fusion plasmas. It fol 
ion viscosities (e.g., gyroviscosities |22l |) must be employed to adequately explain experimental 
observations. The consideration of such closures in our theoretical framework will be the subject of 
future work. 
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